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i .. 

K^ ■ We consider polymers made of magnetic monomers (Ising or Heisenberg-like) 
!>■ 

^^ ■ in a good solvent. These polymers are modeled as self-avoiding walks on a 

(N ' 

(^ ■ cubic lattice, and the ferromagnetic interaction between the spins carried by 

^^ ■ the monomers is short-ranged in space. At low temperature, these polymers 



undergo a magnetic induced first order collapse transition, that we study at 



'^ ■ the mean field level. Contrasting with an ordinary point, there is a strong 

O ■ jump in the polymer density, as well as in its magnetization. In the presence of 



a magnetic field, the collapse temperature increases, while the discontinuities 
decrease. Beyond a multicritical point, the transition becomes second order 
and 0-like. Monte Carlo simulations for the Ising case are in qualitative 
agreement with these results. 
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I. INTRODUCTION 

Under various solvent conditions, a polymer chain can be either swollen or collapsed. In a 
bad solvent, a phase transition between these two states can occur as temperature is varieda. 
This is the so-called O point, the tricritical nature of which has been demonstrated by de 
Gennesa. The effective monomer-monomer attraction results from tracing out the solvent 
degrees of freedom. At the O point, the second virial coefficient of the polymer vanishes. 

In this paper, we consider a different mechanism that also yields attractive monomer- 
monomer interactions. The model we study consists of a polymer chain, where each monomer 
carries a spin S in an external magnetic field ; these spins interact with each other with 
a short-ranged interaction. To be specific, we will consider a self-avoiding walk (SAW) of 
length N, on ad dimensional cubic lattice, with a nearest-neighbor ferromagnetic interaction 
(on the lattice) between the spins of the monomers. 

Several models involving both polymeric and magnetic-like degrees of freedom have been 
introduced in very different contexts. A similar, but somewhat more complicated model, 
was studied to describe secondary structure formation in proteins^: there, the dominant 
interactions between the monomers are of (electric) dipolar nature. A Potts model on a 
SAW was studied as a description of vulcanization inQ; a similar model with quenched 
disorder was studied in the context of secondary structure formation in proteinsH. 

Ising models have also been studied on a fixed SAW geometry, yielding results quite different 
from those presented hereoQ. From an experimental point of view, organic polymeric magnets 
have recently become of interestffl. However, there has been very little interest in their 
conformational changes. 

In the following, we will specialize to the Ising case, and quote some results for the Heisenberg 
case. The partition function of the system reads: 

z=y: y. ^^pi^J:s^^nr,s,+^hJ:s^ (i) 

SAWSi=±l \ ^ ijtj i I 



where J is the exchange energy, /3 = ^ the inverse temperature, and h the external magnetic 
field. The spatial position of monomer i with spin Si is r^. The symbol A^^.' is 1 if {r^r'} 
are nearest-neighbor on the lattice and otherwise. The sums run over all possible SAW 
and all spin configurations. 

This model will be studied along three lines. In section 0, we derive upper bounds for 
the free energy of the model (in zero magnetic field) which suggest a first order transition 
between a swollen paramagnetic and a collapsed magnetized phase. In section |T|, we derive 
a mean-field theory for the general model. We show that indeed for low fields, there is a 
such a line of first-order transitions. At higher fields, this transition becomes continuous. 
The two regimes are separated by a multicritical G point. At this special point, both the 



second and the third virial coefficient vanish. In section |I^, these predictions are tested 
against Monte Carlo simulations in rf = 3 dimensions. The numerical results are consistent 
with the theoretical phase diagram; at low field (where the transition is strongly first order), 
the agreement is even quantitative. 

II. FREE ENERGY BOUNDS 

The physics of the model can be described according to the following simple picture: 

1. At high temperature, since entropy dominates, the chain is swollen. As a result, the 
number of nearest-neighbor contacts is small, and from a magnetic point of view, the 
system is equivalent to a one dimensional Ising model. This simple picture can be 
expressed through the following inequality: 

Z > ZsAW Zi{h) (2) 

where ^saw is the total number of SAW, and Zi{h)is the partition function of the one 
dimensional Ising model in a field h. Using well known results"'Eil: 



^SAW 



/i^ATT-i (3) 



where 7 is a critical exponent, we obtain: 

— <-T\ogfi-T log(e^'^ cosh ph + ^j e'^^-^ sinh^ (3h + e-^^J) (4) 

The best estimate in ci = 3 is // ^ 4.68 for the cubic lattic 

2. At low temperature, the magnetic energy is larger than the entropy loss due to confine- 
ment, and thus the chain collapses. The simplest picture is that of a totally magnetized 
and fully compact system, resulting in the following bound: 

Z > exp{N(3Jd + N(3h) Zhp (5) 

where Z-^p is the entropy of Hamiltonian paths (HP) on the lattice, i.e. fully compact 
SAW on the lattice []. Using a virtually exact upper boundll^ to Zhp, we obtain: 

L<-T\og^-Jd-h (6) 

Note that in two dimensions and in zero magnetic field, it is possible to write a more 
accurate bound by using the exact expression Z2 for the Onsager partition function of 
the Ising model: 

Z>Z2 Znp (7) 

The free energy bounds of equations (^) and (^ are shown as functions of temperature in 
Figure 1, for d = 3 and h = 0. The true free energy lies below these two curves, and their 
intersection is an indication for a possible first order transition between the swollen and 
collapsed phases. 



^A HP is a SAW which passes through each point of the lattice exactly once 



III. MEAN-FIELD THEORY 



By using a Gaussian transform, it is possible to write a field-theoretical representation for 
the model (0): 



Z = 2^ l[[d(pr exp -— - ^ iprA^^l,(pr' + log ^ n cosh{ipr, + (3h) 

r \ '^P'^ {ry} SAW{n} J=l / 



(8) 



As usually, the mean-field theory can be obtained by performing a saddle-point approxima- 
tion on equation (|^). We assume that the chain is confined in a volume V with a monomer 
density p = y. Assuming a translationally invariant field ip, the mean field free energy per 
monomer is 

/ = ^ = -riog2 + ^y.2 _ TlogZsAw - Tlogcosh(y. + ph) (9) 

where q = 2d is the coordination number of the cubic lattice and ^saw is the total number 
of SAW of A^ monomers confined in a volume V. Following ref.tj, it is easily seen that: 

ZsAw - (^) exp {-V{1- p) log(l - p)) (10) 

so that 

/ = -T log 2 + -^<^' - T log ^ + T ^—^ log(l -p)-T log cosh((^ + ph) (11) 

zpjq e p 

This free energy is to be minimized with respect to the field ip and to the volume V occupied 
by the chain, or equivalently to the monomer concentration p. The mean field equations 
read: 

2 

T^q = -P- l°g(l - P) (12) 

(f = (3Jpq tauh{ip + (3h) (13) 

Note that this free energy holds also for a melt of chains, where p is the total monomer 
concentration. 



This set of coupled equations has a high temperature solution p = 0,ip = 0, which describes 
the swollen phase with no magnetization and vanishing monomer concentration, and a low 
temperature solution, which describes a collapsed phase with a finite monomer concentration 
and magnetization. More precisely, for magnetic fields h < ho, there is a first order transition 
(as a function of temperature) between a swollen and a collapsed phase. At higher fields 
h > ho, the transition becomes second order (in fact tricritical). For infinite fields, the 
magnetization saturates, and the model becomes equivalent to the ordinary B point as 
studied in many paperslijii^. 

Expanding (0) and (ITsD , one obtains the equation for the second order line: 




tanh/3/i=W — (14) 



Close to this critical line, the concentration varies as: 

1 I- pjq tanh^ Ph 



(15) 



"^ 2(/5Jg)(l-tanh2/5/i)-| 
and the magnetization per spin is given by 

M ~ tanh/5/i (16) 

and remains finite, whereas the magnetization per unit volume, given by 

m ~ ptajih f3h (17) 

vanishes. 

The phase diagram is shown in Figure 2, with values corresponding to dimension d = 3. 

The first order and continuous transitions are separated by a multicritical point, denoted by 
G on Figure 2. The corresponding temperature and field are Tq = 4.5J and he = 5.926J. 
At zero magnetic field (point A), the transition temperature is Tc = 1.886 J, the critical 
concentration is pc = 0.87 and the critical magnetization per unit volume is rric = 0.87. 



Note that the magnetic susceptibihty x = dM/dh remains finite along the second order 
critical line. 

The same phase diagram holds for the Heisenberg ferromagnet, with a multicritical G point 
at Tg = 2.31J, ho = 5.91J. The zero-field point A is at T, = 0.844J, p, = 0.902 and 
rur = 0.73. 



Using ([13|) to eliminate y9 as a function of p, it is possible to express the free energy (|Tl]) 
only as a function of the monomer concentration p. This yields the virial expansion of the 
free energy. The second virial coefficient vanishes along the second order line (0) , implying 
that the transition is 0-like (i.e. tricritical). At the multicritical G point, both the second 
and third virial coefficients vanish, while the fourth order coefficient is positive. 

IV. MONTE CARLO SIMULATIONS 

The Monte Carlo methodli3 used to compute thermodynamic as well as geometric properties 
of the magnetic chain relies on the multiple Markov chain sampling. A detailed description 
of this method can be found inll3ll^. The implementation we consider for a magnetic chain 
on a three dimensional cubic lattice, can be summarized as follows. 

We start from a single Markov chain at fixed temperature T. The probability n-z){T) of a 
(magnetic) chain configuration V, is given by the Boltzmann distribution txv{T) ~ e^^'yT^)/'^ 
with 

HiV) = -^ E SrKr.S, -hY^S, (18) 

where the thermodynamic variables Si and r, are assigned their D-dependent values. This 
Markov chain is generated by a Metropolis heath bath sampling based on a hybrid algorithm 
for chains with pivotEj a well as local movesc^. Pivot moves are nonlocal moves that assure 
the ergodicity of the algorithm; they operate well in the swollen phase but their efficiency 
deteriorate close to the compact phase. In this respect local moves become essential to speed 



up the converge of the Markov chain. Finally, in addition to the moves that deform the chain, 
an algorithm based on Glauber dynamics is considered to update the spin configuration along 
the chain. For a single Markov chain we typically consider ~ 10^ pivot moves intercalated 
by N/4: local moves and A^ spin updates. 

The multiple Markov chain algorithm is then implemented on the hybrid algorithm described 
above. The idea is to run in parallel a number p (in this work p = 20—25) of Markov chains at 
different temperatures Ti > T2 > ■ ■ ■ > Tp. In practice, this set of temperatures is such that 
the configurations at Tj and at Tj+i have considerable overlap (implying that Tj and T^+i 
are close enough). We let the Markov chains interact by possibly exchanging configurations 
as follows. Two neighboring Markov chains (i.e. with temperatures Tj and Tj+i) are selected 
at random with uniform probability. A trial move is an attempt to swap the two current 
configurations of these Markov chains. If we denore by tik{T) is the Boltzmann probability 
of getting configuration K at temperature T, and Vj and I?j+i are the current states in the 
jth and (j + l)th Markov chain, then we accept the trial move ( i.e. swap T>j and 'C'j+i) 
with probability 

r{V„ Vj+i) = mm 1, \ \ 19 

The whole process is itself a (composite) Markov chain that is ergodic, since the underlying 
Markov chains are themselves ergodic. It turns out that the swapping procedure dramatically 
decreases the correlation times within each Markov chain with little cost in CPU time since, 
in any case, one is interested in obtaining data at many temperaturescB'ta. 

For each multiple Markov chain run we compute estimates, at a discrete set of temperatues 
T, of quantities such as (i) the average energy {£) and specific heat C = - — ^2 of the 
chain. The per monomer quantities will be denoted respectively by E and C (ii) the average 
magnetization per monomer M = j^ Y.i{Si) (iii) the susceptibility x = ^- ^^ addition, as a 
geometric quantity, we consider the mean squared radius of gyration (i?^) of the chain. From 
now on, we will set J = 1, which amounts to give the values of the field and temperature in 
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units of J . 

We have done preliminary simulations at high (T = 10) and low (T = 1) temperature. Our 
results (Figure 3) show that the chain undergoes a swollen to collapsed phase transition, 
in broad agreement with the mean field picture; moreover, the radius of gyration is found 
to vary very little with the magnetic field (note that mean field theory yields pc = 0.87 at 
Tc in zero field). Since a full exploration of the {h,T) plane is difficult, we have restricted 
this paper (i) to a detailed study of the h = transition (ii) to a qualitative study of some 
non zero magnetic field transitions. As mentionned above, the infinite field case corresponds 
to the usual G situation: equation (0) then yields Tg = q = 6, whereas the experimental 
(cubic lattice) value is Tq ~ 3.7t3. As expected in the presence of fiuctuations, the mean 
field parameters (including the location of point G) are not reliable. For small h, the first 
order caracter of the mean field transition will be seen to improve the situation. 

A. Results for h = 0: evidence for a first order transition 

Our results for the specific heat per monomer C and the susceptibility x ci^re respectively 
given in Figures 4 and 7. The spiky character of both contrasts with the rounded specific 
heat of a usual B point. Indeed, finite size scaling theorynll predicts that the peak Cmax of 
the specific heat behaves, in the critical region (i.e. for large enough A^), as 

Cma. ~ N^ (20) 

where a is the critical exponent associated with the temperature divergence of the specific 
heat. Accordingly, the critical temperature shifts from its A^ = cxd value by an amount AT 
given by 

AT ~ N'^ (21) 

At the point, one has a = 0, implying a slow (logarithmic) A^ dependence of Cmax and 
^T ~ j^ (up to a logarithmic factor). On the contrary, a thermal first order transition 
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corresponds to the value a = 1, yielding Cmax ~ N with a much weaker temperature shift 
AT ~ ^. These scaling predictions are to be compared with the results of Figures 5 and 
12. The agreement is satisfactory, even though it is not clear that the largest A^ value, viz. 
A^ = 400, is already in the scaling region. Further evidence for a discontinuous zero field 
transition comes from Figures 6 and 8, where we show the thermal evolution of the average 
magnetization per monomer M and of the radius of gyration. All these results are consistent 
with a first order transition at a critical temperature T^ — 1.80 ± .04, close indeed to the 
mean field value T*^^ ~ 1.88. 

To study the phase coexistence implied by such a transition, we have studied the probability 
distributions of the magnetization M and internal energy E close to the phase transition. 
Figures 9 and 10, obtained for A^ = 300, suggest that the critical distributions P{M) and 
P{E) are fiat, in marked contrast with the usual two peak structure at T^. This two peak 
structure results from the spatial coexistence of the (bulk) phases along a domain wall (more 
generally a (d — 1) interface). In the present case, we have coexistence between phases of 
different dimensionalities, namely a paramagnetic swollen phase and a magnetized collapsed 
phase. It is then clear that the "interface" can be reduced to a point, yielding a "surface" 
tension of order one. This in turn explains the fiat critical distributions of Figures 9 and 10. 
Below the transition, it is interesting to note that the magnetization quickly saturates: a 
closer look at compact chain magnetic conformations shows that the minority domains are 
located on the surface of the globule, and become less and less relevant as A^ grows (Figure 
11). 

B. Tentative studies of the continuous transition in a field 

As previously mentionned, the multicritical point G will be pushed downwards from its mean 
field location. Since the computer search for a precise determination of this point is very 
time consuming, we have adopted the following strategy. We have performed simulations 
for small {h = 0.5 and h = 1) and large {h = 5 and h = 10) magnetic field. 
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The first evidence for a second order transition in large fields conies from Figure 12, where 
the specific heat maximum for h = 10 behaves very differently from its small field values: 
for h = 0.5 and h = 1, one apparently gets the same behaviour as with h = 0, namely 
Cmax ~ N. The existence of a second order transition for h large is corroborated by Figures 
14 and 15. For h = 5, the critical probability distribution is very different from its /i = 
counterpart (Figure 10). A finite size scaling analysis of the data for the same value of the 
field yields, in a rather convincing manner, a second order B like transition at Tc ~ 3.4 
(remember that \im.h->oo Tc{h) ^ 3.7). We therefore obtain he < 5. To get a better estimate 
of he, we have computed the probability distribution P{E) of the internal energy for h = 0.5 
(Figure 13). It clearly interpolates between Figures 10 and 14, but it is not easy to interpret 
the data as representative of a continuous or discontinuous transition. To summarize, we 
have presented evidence for a continuous transition for large h. The precise position of the 
point G is left for future work. 

V. CONCLUSION 

We have seen that ferromagnetic interactions may drive the collapse of a polymer, even in 
a good solvent. This collapse is very sensitive to the presence of an external magnetic field. 
It might be possible to design new polymeric magnetic materials, for which the collapse 
transition is triggered by a magnetic field, at room temperature. 

We have also done preliminary simulations on the two dimensional case (Ising polymer on 
a square lattice): for h = 0, we get a quite abrupt transition around T^ ~ 1.18 (which can 
be compared to the value Tg ^ 1.5c3). Since the critical dimensions associated with the 
point {ip^ theory) and the multicritical G point {ip^ theory) are respectively de = 3 and 
do = ^, one expects fluctuations to be important. Further work is needed to elucidate their 
influence on the mean field phase diagram. 
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Finally, the present model can be generalized to include 

(i) longer range or competing interactions (e.g. ANNNI models). 

(ii) non Ising local variables {0{n) spins, quadrupoles,...). 

(iii) disorder, either in an annealed (BEG-likeEfl) or in a quenched wa;''^ 
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FIG. 1. The high temperature (soHd Hne) and low temperature (dashed hue) free energy bounds 
for d = 3 and h = 0. 
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FIG. 2. Mean- field phase diagram in d = 3 for Ising spins. The dashed hne corresponds to a 
first order transition, and the solid line to a second order (0-like) transition. 
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FIG. 3. Log-Log plot of < i?2 > ^s N, for N = 80, 100, 150, 200, 300, 400. The empty 
and filled symbols correspond respectively to T = 10 and T = 1. Values of the magnetic field are 
h = (o), 0.5 (o), 1 (D) and h = 10 (A). The lines have respective slope 2z/ = 1.194 ± .005 and 
2u = 0.63 ± .04. 
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FIG. 4. Specific heat per monomer vs temperature for h = and iV = 80 (o), 100 (•), 150 (□), 
200 (■), 300 (A), 400 (a). Note the increase of the peak as well as its shape, when A^ increases. 
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FIG. 5. Temperature location ofthe specific heat peak vsl/iV, for Af = 100, 150, 200, 300, 400. 
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FIG. 6. Magnetization per spin vs temperature for h = and A^ = 80 (o), 100 (•), 150 (□), 
200 (■), 300 (A), 400 (a). 
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FIG. 7. Susceptibility vs temperature for /i = and iV = 80 (o), 100 (•), 150 (□), 200 
300 (A), 400 (a). 
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FIG. 8. Squared radius of gyration of the polymer as a function of temperature for h = 0, and 
N = 80 (o), 100 (•), 150 (□), 200 (■), 300 (A), 400 (a). 
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FIG. 9. Probability distribution P{M) of the magnetization per monomer, for three different 
temperatures at /i = and A^ = 300: T = 1.67 < Tc{N) dotted line, T = 1.81 > Tc{N) solid line, 
T = 1.78 ~ Tc{N) thick line. 
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FIG. 10. Probability distribution P{E) of the internal energy per monomer for three different 
temperatures at /i = and A^ = 300: T = 1.67 < Tc{N) dotted line, T = 1.81 > Tc{N) soUd line, 
and T = 1.78 ~ Tc{N) thick line. 
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FIG. 11. Typical chain configuration just below the transition (T = 1.65), for h = 0. Red 
(resp. blue) monomers have up (resp. down) spins. Note that the down spins are on the surface 
of the globule. 
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FIG. 12. Peak of the specific heat per monomer vs N for different values of h: 
h = (o), 0.5 (•), 1 (□), 10 (■). For h = 10, the value of Cmax is clearly not proportional 
toiV. 
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FIG. 13. Probability distribution P{E) of the internal energy per monomer for three different 
temperatures sX h = 0.5 and N = 300: T = 1.74 < Tc{N) dotted line, T = 1.89 > Tc{N) soUd line, 
T = 1.81 ~ Tc{N) thick line. 
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FIG. 14. Probability distribution P{E) of the internal energy per monomer for three different 
temperatures at /i = 5 and iV = 300: T = 1.67 < Tc(iV) dotted line, T = 4 > Tc{N) solid hue, 
T = 2.43 ~ rc(iV) soHd thicker hue. 
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FIG. 15. Search for a second order like transition for h = 5: scaled radius of gyration {-jj-) 
vs temperature for N = 100 (•), 150 (D), 200 (♦), 250 (A), 300 (+), 350 (*), 400 (▼). A crossing 
occurs for T ~ 3.4. 
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